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Abstract. Numerous theoretical and experimental studies have investigated the 
dynamics of cold atoms subjected to time periodic fields. Novel effects dependent on the 
amplitude and frequency of the driving field, such as Coherent Destruction of Tunneling 
have been identified and observed. However, in the last year or so, three distinct types 
of experiments have demonstrated for the first time, interesting behaviour associated 
with the the driving phase: i.e. for systems experiencing a driving field of general form 
V{x)sm{uJt + <fi), different types of large scale oscillations and directed motion were 
observed. We investigate and explain the phenomenon of Super-Bloch Oscillations 
(SBOs) in relation to the other experiments and address the role of initial phase in 
general. We analyse and compare the role of </> in systems with homogeneous forces 
(V'(x) = const), such as cold atoms in shaken or amplitude- modulated optical lattices, 
as well as non-homogeneous forces (V'(x) ^ const), such as the sloshing of atoms in 
driven traps, and clarify the physical origin of the different (^-dependent effects. 
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1. Introduction 

Experiments on cold atoms in optical lattices subject to time-periodic perturbations 
[U 121 |3l IH [5] have provided exceptionally clean demonstrations of coherent phenomena 
predicted in theoretical studies O El [101 HH H21 [13l [14]. Of especial note are 
the phenomena termed Coherent Destruction of Tunnelling (CDT) [TJ El E] or Dynamic 
Localization (DL) [6]. 

In particular, if an oscillating potential V(x, t) = Fx cos ut is applied, it was found 
that the tunnelling amplitudes J of the driven atoms take an effective, renormalized, 
value: 



where d is the lattice constant and j7o denotes a zero-th order ordinary Bessel function. 
In the CDT (u ^> J) regime, there is a complete cessation of tunnelling at the zeros of 
the Bessel function; in the DL regime, the atoms make a periodic excursion but return 
to their original position if jTo ~ [T5] . 

The above showed that even the one-particle dynamics can exhibit non-trivial effects 
arising from matter-wave coherence. The potentials were realized with ultracold atoms 
in shaken optical lattices; Eq.fll]) was demonstrated and investigated experimentally in 
a series of ground-breaking experiments [H [3] . The effects of interactions have also been 
investigated [HI [161 El EIEl EE] such as AC-control of the Mott-Insulator to Superffuid 
phase transition, including even triangular optical lattices [20] . 

The experiments above identified effects which depend on the frequency and 
amplitude of the driving field, not its initial phase. However, in the last year or so, 
three distinct experiments on periodically driven cold atoms have identified effects due 
to the phase; in other words, for driving of form V(x,t) oc sm(ut + <p) , the dynamics was 
found to depend strongly on </>. For example, for atoms in shaken lattices, [211 [22] were 
able to realize both directed motion as well as large oscillations, occurring over hundreds 
of sites, which were termed "Super-Bloch Oscillations" (SBOs) [22J. These SBOs were 
analyzed in |23j, including also the effects of weak interactions (mean- field regime). The 
general assumption in theoretical analyses has been that the external field-dependence 
of the wavepacket mean group velocity in these cases is entirely contained in the Bessel 
function argument, with no dependence on 0. However, the experiment indicates that 
there is a strong dependence on this phase. 

In a second example, [21], the tunneling amplitude of the lattice was modulated. 
The result was a global motion (non-zero velocity of the atomic wavepacket), strongly 
dependent on <fi. The third example did not involve optical lattices or tunneling 
but looked at a BEC "sloshing" in a time-averaged orbiting potential (TOP) trap. 
Experiments showed [22] using a 2-dimensional set up, that the amplitude of oscillations 
of atoms in traps has a strong dependence on phase. Related theoretical work on the 
atoms oscillating in traps was carried out by J2SI [2Z] ■ 
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The atoms in traps are also driven with potentials of the form (in 1-dimension) 
V(x) sin (cut +0). Yet they differ in essential ways. Because the potential is not linear, the 
force is not uniform; the time-averaged motion becomes strongly coupled to the phase; 
the trap analysis always requires application of an adiabatic separation between fast u> 
and slow mean motion. In addition, they also differ from the optical lattice systems 
because the latter are band-Hamiltonians, with an effective kinetic energy —Jcosp (in 
a classical- image Hamiltonian) . Such differences mean that effects due to <ft are more 
subtle in the lattices; thus for the latter, phase effects have been neglected for several 
years, while they quickly became apparent in the trap systems. 

However, at present, the physical role of the initial phase in all these experiments 
does not seem too clear. Recent analyses consider the effect of the phase in terms of 
a shift in the origin of time for shaken lattices [4] or that the effects arise due to the 
abrupt amplitude jump which occurs when the field is initially switched-on [25] at t ~ 0. 

We have recently found [28] , that in order to explain the experiments [211 122] , a 
phase correction $(-F) = F w cos 0— n(0+|) must be considered in an effective dispersion 
relation. This explains a field dependent shift in the phase of SBOs. It accounts also for 
a separate regime of directed motion. Here we analyse these SBOs in relation to other 
comparable experiments with cold atoms which have also observed other types of global 
motion and large scale oscillations such as trap oscillations or amplitude-modulated 
dynamics. More broadly, the aim of this work is to clarify the role of the initial phase in 
such systems and to consider the validity of proposed models and assumptions in these 
related experiments. 

It is worth stressing that for shaken lattices and SBOs, the strongest effects occur 
when the system evolves from p = with a smooth pure sin ut drive, where there is no 
amplitude jump or even a phase-jump at t — 0. Yet below we consider only oscillations 
where the phase is well-defined in each oscillation. Thus <fr has to be set on a timescale 
fast compared with T = 2tt/u. Variants with a slower ramping- up procedure depend 
on the particular protocol and we do not consider these. 

Below we consider systems with homogeneous forces (V'(x) = const) separately 
from those for which the forces are position dependent. 

2. Inhomogeneous systems: Atoms in Traps 

We consider first the classical dynamics for a Hamiltonian of form: 



Hq(p) = — J cosp (for a one-body image Hamiltonian for a particle in a band). One may 
also consider an additional time independent potential, i.e. H (p,x) = H (p) + Vq(x) 
but for the essential physics, this is not necessary for now. 
Thus there is a classical force: 




p= -V\x)F(t). 



(3) 
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If one can neglect the time evolution of x(t) in the interval [0 : t], we have a momentum 
shift on the initial momentum p : 

p(t)=p -V'(x ) fF{t)dt, (4) 

Jo 

where Xq = x(t = 0) 

One instance is the situation, well-studied experimentally [291 Ell EH E2} |3U E] , 
of ^-kicked cold atoms where an optical lattice potential is abruptly switched-on for a 
very short time and where F(t) = S n 5(t — nT). For a single kick, and T — 1, the 
shift is p = po — V'{x). There is an instantaneous impulse applied to the atoms. It 
is well-known that, provided the switch on is sufficiently fast and the pulse is of short 
duration, the dynamics is insensitive to the pulse shape. However, the shift is position 
dependent; for repeated kicks, this can have drastic consequences including chaos, for 
which kicked atoms leading paradigm. 

Another instance of a position-dependent driving, corresponds to atoms in traps 
[26| [27] . where V(x) represents, for example, a harmonic trap of frequency Q and 
F(t) = sin(a;t + 0) represents a sinusoidal drive. In the case u ^> Q one can consider 
a period-averaged mean motion for the slow oscillation in the well, with coordinate 
X(t), as well as a "micromotion" characterized by the faster u timescale. Adiabatic 
separability means that one may consider these motions separately. In [26j [27] an 
effective shift of the mean- momentum P Q = mX(0) was identified. In the present 
notation, one would obtain an effective phase-dependent shift of the period-averaged 
momentum: 

pi =Po + Y!M COS(j) . (5) 

A mean-momentum shift 

AP(x ) = ^^cos0 (6) 

CO 

is zero for = 7r/2 or = 3n/2. But it is maximal for = 0,7r the case where 
the F(t) = ± sin ut; for this case, the driving field starts smoothly from zero. This 
clearly indicates that this case represents essentially different physics from the ^-kicked 
particles. The shift is unrelated to the abruptness of the change in the Hamiltonian 
amplitude at t = 0, but rather to a "slippage" discussed in [26], [27] between the slow 
mean motion and the fast "micromotion" . It makes little sense to define period-averaged 
motion on a timescale shorter than of order T = 2tt/uj: the system requires a finite time 
to "notice" that the mean momentum has shifted; for the 5-kicked particle, in contrast, 
there is a real impulse, applied instantaneously on a timescale <C T. A formal way of 
stating this is that a 8(t) function implies a finite impulse at t ~ 0, while the Heaviside 
step function switch-on of a sin(a;£ + 0) drive gives zero impulse at t ~ and is thus 
not significant to the dynamics. 

Nonetheless, the shift AP(io) has real physical implications as it can reduce or 
increase the kinetic energy of the mean motion, depending on its sign. It leads to an 
increase or decrease of the amplitude of the harmonic oscillations of the particle in the 
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trap (averaged over one period). Ridinger and Davidson [2S1[27], for instance considered 
these and other examples of traps (Eqs. (15) and (16) therein). 

These effects were demonstrated experimentally in atoms sloshing in a TOP trap 
|25j . The experiments however, used a 2-dimensional set-up. Even for = 0, a rotating 
field (B cos cut, B sin cot) is applied in the (x,y) directions. There is never a situation 
where a pure sine drive is applied and there is always an initial "jump" in field amplitude 
in at least one coordinate. Thus a convenient view (Ref. [25J) is to reset the large scale 
oscillation (macromotion) amplitude coordinates at the end of the (near-instantaneous) 
switch-on of the field. The macromotion is thus attributed to the sudden switch-on. 
However, one should stress that the underlying physics can emerges only on the timescale 
of T, rather than being due to the abrupt initial switch-on procedure. 



3. Homogeneous systems: shaken optical lattices 

The shaken optical lattices represent a different class of Hamiltonians. The 
corresponding classical Hamiltonian is: 

H(x, t) = — Jcosp — Fxsm(cot + 0). (7) 

A key difference is that V'(x) = F is independent of position. While the particle 
does display Bloch Oscillations (SBOs) with the period Tb, there is no separation of 
timescales since Tg ~ T. There is no restriction to high frequency driving. Further, 
there is a decoupling between the instantaneous phase of the Bloch Oscillations and the 
phase, since the initial position x(t — 0) is immaterial: the force is uniform. The initial 
wavepacket is even delocalized over several lattice sites. From the classical equations 
of motion, p = p + F J * sin (cut' + (f))dt' and a effective momentum shift arises from the 
integrand at t = 0: 

F F 

pit) = po H cos0 cos(cut + 0). (8) 

CO CO 

The effective shift is, however, position independent: p' Q = p + ^cos0. It is possible 
to then consider the effect of the phase as simply a corresponding displacement in the 
origin of time as suggested in jl], but only in the momentum shifted frame of p' . 

Once again, for the completely smooth switching on of the driving Hamiltonian 
(pure sin cut drive) the shift is a maximum. And here it is not useful to think of an 
initial "jump" either in amplitude or phase. Taking the experimental value po = and 
= 0, 

F 

(p(t)) = — (1 - cos cot). (9) 

CO 

We see that the momentum now oscillates about a shifted non-zero value, but the 
process is not related to the abruptness of the switch-on, which is perfectly smooth. The 
shift in the average momentum becomes apparent only on the timescale of T. This is 
very different from the 5-kicked particle, where a real instantaneous impulse is imparted. 
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It is also different from the trap, where through V'(xq), the position at the instant the 
phase is set is important. 

One may solve Hamilton's classical equations of motion to obtain a period averaged 
group velocity: 

v g ~ -JJoiF^) sin (p + cos 0) , (10) 
where F w = F/ui. 

In the shaken lattice, there is no dependence on initial position xq. Thus we 
conclude that the directed motion is not related to the sudden switch-on of the field. It 
is best understood as a process more akin to ratchet physics and the de-symmetrization 
of the drive. The time-asymmetric sin cut drive generates a ratchet current (since the 
average force is zero, there is no net bias: thus this fulfils one of the criteria for ratchet 
motion, though not the stricter one of "rectification of fluctuations" , so should perhaps 
be termed simply "directed motion"). The symmetric drive cos cut produces no ratchet 
motion, while intermediate values of produce intermediate behaviour. The final 
current is similar to that obtained in the linear-ramping study of Ref . [35] . 

The above discussion explains only directed motion and its dependence. For 
Super-Bloch Oscillations (SBOs), one has an additional static field. One must consider 
in addition the quantum resonance between u and the gap between energy levels 
in adjacent wells. We thus consider a quantum treatment, within a Floquet-theory 
framework (full details are in [28]). 

We consider the total Hamiltonian: H(t) = H + Hp(t), where 

^o = -^E( c L c ^ + H - c -)- (ii) 

H corresponds to the non-interacting limit of a variety of Hamiltonians with nearest- 
neighbour hopping (Hubbard, Bose-Hubbard, magnons in Heisenberg spin chains, etc.) 
It represents any spatially periodic potential characterized by energy eigenf unctions 
4> m k, with band index m and wavenumber k, thus H (f) m k = E m (k)4> m k- We restrict our 
one-particle problem to the lowest band m — 1; taking E(k) = E m= i(k), the energy 
dispersion E(k) = — J cos kd. 

We take Hp[t) = —f(t)x where fit) = F + F sm(ut + 0); in general it comprises 
both a static field and a sinusoidally oscillating field with an arbitrary phase 0. The 
result of the driving is a time-dependent wavenumber [I] : 

1 /"* F F 
q k (t) = k + - drf(r) = k + —cos cos(wt + 0) (12) 

h Jo <jj <jj 

The stationary states of the system are its Floquet states, the analogues of 
Bloch waves in a temporally periodic system. They are given by ipk{x,t) = 
Uk(x,t) exp —^e(k)t , where Uk(x,t) = Uk(x,t + T), and the period T = 2h/oj. The 
non-periodic phase-term is characterized by the so-called "quasienergy" e(k). 

The presence of the static linear field in general destroys the band dynamics; 
however, here we consider the so-called resonant driving case, for which F^d = nhu, 
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(where the driving compensates for the energy offset between neighbouring wells in the 
lattice, restoring the band structure). In this case, other studies, while neglecting </>, 
have found that the renormalization of tunnelling is n— dependent: 

J cS oc JJ n (^j , (13) 

where d is the lattice constant and J n denotes an ordinary Bessel function and 
n — 0, 1, 2... Below we take % — 1 and F u = Fd/ui. 

For the purposes of calculating to group velocities, we define an effective quasienergy 
dispersion: 

t(k) = — F cos{q k (t)d)dt (14) 
1 Jo 

by means of a period-average. 

Below we also consider the case of slight detuning for which F Q d = (n + S)u, with 
5 <C 1, associated with SBOs; however, for the slight-detuning case 5 ^ 0, we assume 
that the time-dependence due to the Sut/d remains negligible over one period T. Thus 
we take it out of the period-average entirely and, as shown in [28], Eq. ( fl4l becomes: 



e(k) ~ -JJniK) cos 



7T 

kd + d~ujt + F^ cos <p — n(<p + — ) 



(15) 



The above represents an effective dispersion relation, but which oscillates slowly in time 
with a period T SB o = 2n/(5u)) ^> T B where T B oc 1/F is the Bloch period. They 
correspond to the "Super-Bloch Oscillations" investigated by Refs. [131 [231 Ell 122] . 
Even at resonance 5u = 0, Eq. ( Ti~5"]) differs from previous expressions by the phase-shifts 
$(F) = F w cos0 — n(<fi + |). Evaluating Eq. (fT5|) for the n = 1 case, we obtain a 
period-averaged group velocity: 
de_ 
dk 

The above expression accounts for a range of experimental features. However, the 
most interesting one is that the SBOs begin with a field-dependent phase. Figure [H 
demonstrates this for experimental field values F^ = 0.08, 0.15, 0.76 and 1.52 (com- 
pared with one-particle Hubbard numerics). The graphs show that Eq. (fl6l) reproduces 
quite well the experimental values of Ref. [22]. They show clearly the displacement 
of the first maximum, seen in the experiment as well as the order of magnitude varia- 
tion in amplitude. Such a field- dependent shift is also apparent in the results of Ref. |21j. 



v g = — ~ — Jd cos (kd + F^ cos 4> — 4> + J\{FJ) ■ (16) 



Finally, we consider the case of experiments with amplitude modulated lattices, 
which do not show renormalisation of tunneling, but can also produce directed motion 
[21] . The experiments also used a static linear field Fqx. An experimental group velocity 
v g oc J cos (kd — 0), for the n = 1 resonance case F d = hu, was measured. Curiously, 
the effective momentum shift F^ cos<p seen in other driving systems is no longer present. 

We find that the observed behavior in [24J is consistent with an effective classical 
Hamiltonian for the homogeneous force system: 

H(x, t) = - J[l + a sm(ut + </>)] cosp + F x. (17) 
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Figure 1. (a) and (b) One-particle solutions of the Hubbard Hamiltonian with 
the driving term for L = 30 lattice sites, showing the field-dependence of the phase 
of Super-Bloch Oscillations [35]. Upper/lower panel correspond to F u = 1.52 and 
0.15; both have (/) = %. For illustration purposes, the theoretical amplitudes were 
equalized by choosing J- values which equalize J c s — JJx{F u ). (c) and (b) Show 
that good agreement is obtained between experiment and the Hubbard numerics, as 
well as the analytical formula [35J: integrating Eq. (fTo]) (dashed lines) reproduces well 
both the amplitude and phase of SBO experimental data (symbols) of Ref. [22], using 
8u) = 2tt/1000 (= -1 Hz detuning), (f> = it and Jd/(h5uj) = 90 fim. 



Solving for Hamilton's classical equations of motion, we see that p(t) = p^—F^t no longer 
has a phase dependence. For simplicity, we take h = d = 1 and thus Fq = u. Then we 
easily see that the period averaged velocity (Ax) /T = ~ Jq x(t)dt = —^Jacos(p — 0). 
In other words, even with zero momentum shift, in this case a non-zero phase correction 
can arise. 

4. Conclusion 

In this study, we have attempted to clarify the role of driving phase in the dynamics of 
homogeneous and inhomogeneous-force systems. To this end we have compared recent 
experimental results on homogeneous systems (shaken lattices) and inhomogeneous 
systems (atoms oscillating in traps). 

It is tempting to draw analogies between the "slow" SBOs with the macromotion or 
"sloshing" seen in traps [25]. However, in the latter, the effect of the phase is manifested 
directly in the amplitude of the oscillations in the trap. In the shaken lattices, the 
manifestation of the phase is more subtle. The amplitude of SBOs is independent of 
0, only their phase is affected. Unlike the trap systems, their time-dependence does 
not arise from the period-average integral. Nonetheless, it seems clear that in order to 
account for the experiments it is necessary to include the phase- correct ion <&(F) in an 
effective dispersion, used to obtain group velocities. 

What all these systems do have in common is that the main effect of the initial 
phase is to modify the global centre of mass dynamics of the wavepacket. This is in 
contrast to previous experiments on, for example, DL, which showed diffusive spreading 
of the wavepacket rather than large-scale motion, on scales of many sites. 
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